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ABSTRACT 

Linear  programming  techniques  are  becoming  of  greater 
importance  because  the  use  of  computerization  has  increased 
the  fields  for  applications  for  linear  programs.   The  primal- 
dual  algorithm,  in  which  the  constraints  are  added  one  at  a 
time,  is  investigated  as  a  possible  faster  solution  method. 
A  computer  program  was  developed  to  compare  this  method  with 
the  standard  primal-dual  algorithm  using  the  full  set  of 
constraints  at  one  time.   Several  random  problems  were  solved 
using  these  two  methods,  and  the  results  indicated  a  signif- 
icant improvement  in  the  solution  time  by  the  use  of  adding 
the  constraints  one  at  a  time. 
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I .   INTRODUCTION 

New  linear  programming  algorithms  have  been  developed  to 
reduce  the  computational  time  in  solving  linear  programs.   The 
purpose  of  this  thesis  is  to  investigate  the  merits  of  one 
such  new  algorithm.   This  method  consists  of  introducing  the 
constraint  equations  one  at  a  time.   After  each  constraint  is 
added,  the  "smaller",  or  submatrix,  problem  is  solved  using 
the  primal-dual  algorithm.   This  continues  until  all  con- 
straints have  been  added  and  a  solution  is  obtained. 

The  rationale  for  this  approach  is  that  small  matrices 
are  used  in  the  initial  stages  of  solving  the  linear  program; 
the  size  of  the  matrices  increases  only  when  additional  con- 
straints are  introduced.   If  the  number  of  iterations  used  in 
this  method  is  not  significantly  different  from  the  number  of 
iterations  used  with  the  full  matrices,  the  manipulation  of 
the  smaller  matrices  in  the  initial  stages  will  reduce  the 
solution  time. 
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n 
A 


-1 


P  . 


p  . 

ai 


B 


s  . 


s  . 

D 


II.   NOTATION 
number  of  constraint  equations, 
number  of  legitimate  variables. 

mxn  matric  of  coefficients  of  the  constraint  equations 
with  elements  a . . ,  i  =  1 , . . . ,m,  j  =  1,. .. ,n  . 
mxm  matrix  of  the  basis  vectors, 
inverse  of  the  basis. 

mxl  column  vector  which  is  the  j    column  of  A. 

a 


i.e..  P  .  = 


ID 


m: 


/  D  —  J-  f  •  •  •  f  n  . 


.th 


mxl  column  vector  associated  with  the  i    artificial 

variable,  i  =  l,...,m  . 

mxl  column  vector  with  elements  x. .,  where  X.  =  P   P.. 

iD         3  1 

nxl  column  vector  with  elements  c  .  which  are  the  costs 
of  the  legitimate  variables. 

mxl  column  vector  with  elements  b.  which  are  the  right- 
hand  sides  of  the  constraint  equations, 
dual  slack  variables. 

dual  slack  variables  after  a  dual  iteration. 

mxl  column  vector  whose  elements  are  the  coefficients 
of  the  basis  variables  of  the  added  constraint 
equation. 


Notation  used  with  tableau: 


p 

B 

PC 

I    Pl    P2   P3   P4   P5 

p  n  p  i 

aO   al 

Pa0 
P3 

bQ-4(  lj  1/2   1/3   0   5/6   1 
4    0    1/2   2/3   1   1/6   0 

1   -1/6 
0    1/6 

z  .-c  . 
]  3 

s  . 

: 

s  . 

3 

4-bQ 

4bo 

16 

-1 
4 
0 

-1/2  -1/3   0  -5/6  -1 
2     4    0    4    4 
0    8/3   0   2/3   0 

bQ  -  4 


=  -4/-1  =  4 


E3 


the  vector  to  be  introduced  into  the  basis,  e.g., 

in  this  tableau  PQ  will  be  introduced. 

pivot  element  for  a  primal  iteration,  i.e.,  the 

6  criterion  is  0  =  min(x._./x,  .)  such  that  x.  .  >  0 

lB'  ii  i] 

l 

pivot  element  for  a  dual  iteration,  i.e.,  the 
9-criterion  is  9  =  min (-s ,/z  ,-c . )  such  that 

j      3   D   : 

z,-c,  <  0 . 
3       3 


III.   FORMULATION  OF  THE  PROBLEM  AND  SOLUTION  PROCEDURE 
The  general  linear  programming  problem  is  to  maximize 

n 

z  =     y     c .x  . 
3-1  :  3 

subject  to 

n 

I      a. .x.  =  b.    0,  i  ■  1,. .. ,m  ,        (1) 
j=l   l:  :     x  ~ 


and 


x.    0 ,  j  =  1 , . . . ,n  . 


The  modified  primal  uses  an  additional  constraint 

n 

x 
3 


xo  +  .1    xi  =  V 


j=l 

where  the  cost  of  xQ  is  zero  and  bQ  is  arbitrarily  large,  so 
that  for  xQ  >  0  the  constraint  adds  no  additional  restriction 
on  (1). 

The  modified  primal  is  written  to  maximize 

n 


j-1 


]  3 


subject  to 


n 


X0  +  I      Xj  ■  b0  ' 

n 

)   a..x.  =  b..  i  =  l....,m  . 
j£l   x^  ^     i 


(2) 


and 


x.  >  0,  j  =  i,...,n  . 
From  (2)  we  can  write  the  modified  dual  with  slack  var- 
iables, s.,  j  =  0,1,..., n,  added.   That  is  to  minimize 

m 

wr.bA  +   7   w.b. 
0  0    .  L ,   11 
i=l 

8 


subject  to 


w0  '  s0  =  °' 

m 

wn+   I   w.a.-s.=c.,j=l,,..,n,       (3) 
0    i^i   1  ij     ]     j   J 

and 

w.  unrestricted  for  i  =  0.1..... m. 
l 

The  starting  feasible  solution  to  the  primal  dual  al- 
gorithm is  w.  =  0,  i  =  l,...,m,  and  wn  =  max  (c . , 0) . 

j     : 

For  an  optimal  solution  the  complementary  slackness  con- 
dition must  hold.   That  is 

n 

s  x  +  I      s.x   =  0  .  (4) 

«  u    j=1  3    3 

Adding  artificial  variables,  x  .,  i  =  0,1,...  ,111,  to  (2) 

ax 

with  the  cost  of  the  artificial  variables  set  to  -1  and  the 
cost  of  the  legitimate  variables  set  to  zero,  the  extended 
primal  can  be  written  as 
maximize 


m 


-  x  n    -      )      x  . 
L ,   ai 


i=l 
subject  to 


n 


x~  +   J   x.+xA     =bA, 


n 


y   a. .x .  +  x 

i=l   ^  ^ 


=  b  ,  ,  i  =  l,...,m, 
ai     1        tii 


and 


x.    0,  j  =  l,...,n,  and  x  ,    0,  i=0,...,m. 
1  —  ai  — 


Solving  the  extended  primal  is  similar  to  using  a  Phase 
I  Revised  Simplex  method.       However,  in  the  primal-dual 
algorithm,  when  Phase  I  ends  the  linear  program  is  solved 
because  complementary  slackness  is  maintained  throughout  the 
solution  procedure. 

If  a  new  constraint  is  added  to  the  tableau,  complemen- 
tary slackness  is  maintained  without  changing  the  dual  slack 
variables.   This  can  be  shown  as  follows: 

Assume  we  have  a  feasible  solution  to  the  problem 
with  k  constraint  equations.   This  means  that  s.  =  0 
for  all  j  such  that  P.  S  P  and  x.  =  0  for  all  j  such 
that  P.  ^  P.   These  conditions  imply  that  complemen- 
tary slackness  is  maintained,  and  that  the  modified 

dual  also  has  a  feasible  solution. 

st 
Now  we  add  the  k  +  1    constraint  which  introduces 

a  new  dual  variable,  w,  , ,  but  no  new  dual  slack 

variables,  s.,  j  =  D,l,...,n.   We  need  a  feasible 

solution  to  the  modified  dual  for  the  enlarged 

system.   Observe  that  we  have  a  feasible  solution  if 

we  set  w,  ,  =  0  since  then  the  s.,  j  =  0,l,...,n, 

remain  unchanged.   In  particular,  s.  =  0  for  all  j 

such  that  P.  6  P,  that  is,  for  the  legitimate  var- 

iables.   Also,  x.  =  0  for  all  j  such  that  P.  ^  P, 

which  implies  that  we  have  maintained  complementary 

slackness . 

It  is  worth  noting  that  the  z.-c.,  j  =  0,1,..., n,  must  be 

recalculated  since  the  new  constraint  which  is  added  to  the 
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extended  primal  starts  with  its  artificial  variable,  x   .  -  , 

Si  i  K."t  _l 

in  the  basis  with  its  cost  set  at  -1. 

P    0  0 

The  new  basis  is  P  =         ,  where  P   .  , ,  = 

_m  a,k+l 

d1   1  1 

is  the  artificial  vector  associated  with  x  ,  , . .   Now  we  can 

a,k+l 

solve  the  new  k  +  1  system  using  the  primal-dual  alqorithm 
since  complementary  slackness  has  been  maintained. 

An  optimal  solution  exists  if  and  only  if  the  following 
criteria  are  satisfied: 

1.  z.-c.  >  0  for  j  =  0,1,. ...n, 

J       1    ~ 

2.  z  -c_.  =  0,  and 

a      a 

3.  xQ  >  0. 

The  solution  procedure  is  as  follows: 

The  first  tableau  is  set  up  using  the  first  two 

constraints  of  the  extended  primal  and  the  starting 

solution  to  the  modified  dual,  which  implies  that 

at  least  one  s.  =  0. 
3 

Step  1.   Is  there  a  j,  say  j~,  such  that  s.Q  =  0  and 

Z30-Cj0  <  0? 

a.  Yes,   Go  to  2. 

b.  No.   Go  to  3. 

Step  2 .   Introduce  P0  into  the  basis  using  the  minimum 

0-criterion  and  a  primal  iteration.   Since  the  extended 

primal  is  bounded  a  pivot  will  always  exist.   Note  that  the 

s.  remain  unchanged  for  all  j.   Go  to  1 . 

Step  3.   Is  z.-c.  <  0  for  some  i? 
1      1  J 

a.  Yes.   Go  to  4. 

b.  No,   Go  to  5. 

11 


Step  4.   Use  the  minimum  6-criterion.   Is  6  bounded? 

a.  Yes.   Perform  a  dual  iteration  to  compute  a  new 

set  of  s.'s,  say  s.'s.   Go  to  1. 
1  x      1 

b.  No.   The  linear  program  has  no  feasible  solution. 
Stop. 

Step  5.  Is  ZB_CB  <  0? 

a.  Yes.   The  linear  program  has  no  feasible  solution. 
Stop. 

b.  No.   Go  to  6. 

Step  6.   Have  all  of  the  constraints  been  added? 

a.  Yes.   Go  to  9. 

b.  No.   Go  to  7. 

st 
Step  7.   Introduce  the  next  restraint,  say  the  k  +  1 

Place  the  artificial  vector  P   .  . ,  in  the  basis.   Compute 

a,k+l 

xB,k+r    Is  xB,k+i  ±  0? 

a.  Yes.   Go  to  8. 

b.  No.   Multiply  all  coefficients  of  the  new  constraint, 
except  for  the  artificial  variable,  by  -1.   This  assures  that 
x   .  .  _>  0  and  the  artificial  variable  is  non-negative.   Go 
to  8. 

Step  8.   For  the  system  with  k  +  1  restraints,  compute  the 
new  values  of  z.-c.  for  j  =  0,1,..., n.   Go  to  1. 
Step  9.   Is  xQ  =  0? 

a.  Yes.   The  linear  program  is  unbounded.   Stop. 

b.  No.   An  optimal  solution  has  been  found.   Stop. 
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IV.   SAMPLE  PROBLEM 
Consider  the  following  example: 
maximize 


z  =  2x,  +  4x~ 


subject  to 


and 


3x,  +  4Xp  +   6x_.  +  x.       =  24  , 


4x,  +  3x„  +  12x..       +  x,-  =  24  , 


x.  +   x_  +   4x~  =   8  , 


x  .  >_  0,  j  =  1,.  ...  ,5. 


Then  the  extended  primal  is 
maximize 


Xa0    xal  "  Xa2    Xa3 


subject  to 


xfi  +  x,  +   x~  +    x^  +  x,    xr    x 


3x,  +  4x~  +  6x~,   +   x. 


4x1  +  3x„  +  12x3 


x,  +   x~  +   4x^ 


,  -5  T 

aO                                              =  bQ, 

+   xal                                =   24, 

♦   x5 

+    xa2                   =    24, 

»    xa3   =      8, 

x.  >  0  for  j  =  0.....5,  and  x    >  0  for  i  =  0,...,3  . 

J     —  J  r  i       r  ai  — 

The  dual  slack  variables  are  sQ  =  max(c. ,0)  =  4  with  j  =  3. 

j    1 

Then  s_,  =  xnbA  =  4bn,  and  s.  =  s-.  -  c   for  j  =  1,...,5,  so 
B     00      0        j     0     -)      J 

that  s,  =2,  s„  =  4,  s -,   =    0  f  s.    =    4,  and  s5  =  4. 
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The  starting  tableau,  using  the  first  original  constraint, 


is 


p 

B 

P0 

Pl 

P2 

P3 

P4 

P5 

Pa0 

Pal 

Pa0 

bo 

1 

1 

1 

1 

1 

1 

1 

0 

Pal 

24 

0 

3 

4 

£ 

)  1 

0 

0 

1 

-bQ-24 

-1 

-4 

-5 

-2 

-1 



z  .-c  . 
3      3 

bzJ 

___ 

s  . 
3 

4bo 

4 

2 

4 

0 

4 

4 





6  =  24/6  =  4 


From  step  1,  we  see  that  s3  =  0  and  z,-c.  <  0.   Using 

the  minimum  6-criterion  (as  discussed  in  section  III)  in  step 

2,  we  introduce  P.,  into  the  basis  and  remove  P  ,  from  the 

•j  El  A. 

basis . 

Since  there  is  no  jn  for  which  s . A  =  0  and  z.n-c.n  <  0, 

J0  jO  ]0   jO 

but  z.-c.  <  0  for  several  values  of  j,  we  arrive  at  step  4. 
D   D 

/\  /\ 

Using  the  minimum  9-criterion  a  new  set  of  s.'s,  called  s.'s 

are  calculated. 


p 

B 

P0   Pl 

P2 

P3 

P4   P5 

Pa0  Pal 

Pa0 

b0-4 

01/2 

1/3 

0 

5/6   1 

1   -1/6 

P3 

4 

0   1/2 

2/3 

1 

1/6   0 

0    1/6 

z  .-c  . 

-bQ+4 

0-1/2 

-1/3 

0 

-5/6  -1 



s  . 

4b0 

4    2 

4 

0 

4    4 



s  . 
3 

16 

0    0 

8/3 

0 

2/3   0 



=  V4 


e=  -4/-1  =  4 


Now  s  =  0  and  zQ-cQ  <  0  so,  from  step  2,  we  introduce 

PA  into  the  basis  and  remove  P  n    from  the  basis, 
u  au 
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p 

B 

P0 

pl 

P2 

P3 

P4 

P5 

PaO 

Pal 

po 

bQ-4 

1 

1/2 

1/3 

0 

5/6 

1 

1 

-1/6 

P3 

4 

0 

1/2 

2/3 

1 

1/6 

0 

0 

1/6 

z  .-c  . 

3       3 

0 

0 

0 

0 

0 

0 

0 





s  . 
3 

16 

0 

0 

8/3 

0 

2/3 

0 





From  the  above  tableau  we  trace  through  steps  lb,  3b, 
5b,  6b,  and  arrive  at  step  7.   Note  that  we  have  obtained  an 
optimal  solution  to  the  subproblem  with  one  constraint.   In 

step  7  we  introduce  the  second  constraint.   Since  P~  and  P- 

-T 
are  basis  vectors,  d   =  (a2Q,a23)  =  (0,12);  the  new  basis 

consists  of  Pfl ,  P3  and  P  2«   With  this  basis  we  find  that 
x  2  =  -24  <  0  so  that  step  7b  must  be  used.   The  second  re- 
straint is  replaced  by 

-41-30-120-xc+x0=  -24  , 
xl     x2      x3     5     a2 

which  is  used  throughout  the  remainder  of  the  solution  proce- 
dure.  Note  that  now  dT  =  (0,-12)  and  xB2  =  24  >  0.   New 
values  of  x.-c.  for  j  =  0,1,..., n  are  computed  (step  8),  and 
the  new  tableau  is 


p 

B 

po 

Pl 

P   P    P   P 

a0    al   a2 

po 

P3 
Pa2 

bQ-4 

4 
24 

1  1/2  1/3   0  5/6   1 
0  a/2)  2/3   1  1/6   0 
0   2    5    0   2-1 

1   -1/6   0 
0    1/6   0 
0     2    1 

z  .-c  . 

:  : 

s  . 
3 

-24 
16 

0 
0 

-2 

0 

-5    0-2    1 
8/3   0  2/3   0 



+ 

1/2 


=  8 
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Since  s,  =  0  and  z,-c,  <  0  (step  1)  we  go  to  step  2. 
Using  a  primal  iteration,  we  introduce  P,  into  the  basis  and 


eliminate  3 

33  f 

:rom 

the 

basis . 

P 

B 

po 

Pl 

P2 

P3 

P4   P5 

P  n    P  1    P  o 

aO   al    a2 

po 

brt-8 

0 

1 

0 

-1/3 

-1 

2/3   1 

1   -1/3   0 

pl 

8 

0 

1 

4/3 

2 

1/3   0 

0    1/3   0 

Pa2 

8 

0 

0 

7/3 

-4 

(4/3)  -1 

0    4/3   1 

-8 
16 

0 
0 

0 
0 

-7/3 
8/3 

4 
0 

z  .-c  . 

-4/3 

1 

s  . 

2/3 

0 



s  , 
D 

12 

0 

0 

3/2 

2 

0   1/2 



4/3 


=  6 


"  Z473  "  1/2 


Now  z.-c.  >  0  for  all  j  for  which  s.  =  0.   From  steps  lb, 

3a,  and  4a,  a  new  set  of  s.'s  are  computed.   Then  s.    becomes 

zero  and  z.-c.  <  0  so ,  from  step  2,  P.  enters  the  basis  and 

P  ~  is  removed  from  the  basis. 
a2 


p 

B 

P0 

pl 

P2 

P3 

P4 

P5 

Pa0 

Pal 

p  0 
a2 

po 

bQ-12 

1 

0 

-3/2 

1 

0 

3/2 

1 

-1 

-1/2 

pl 

6 

0 

1 

3/4 

3 

0 

1/4 

0 

0 

-1/4 

P4 

6 

0 

0 

7/4 

-3 

1 

-3/4 

0 

1 

3/4 

z.-c, 

J     J 

0 

0 

0 

0 

0 

0 

0 







s  . 

12 

0 

0 

3/2 

2 

0 

1/2 







Optimality  has  now  been  obtained  with  the  second  restraint 
added.   Following  steps  lb,  3b,  5b,  6b,  and  7,  we  introduce 
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the  third  and  final  restraint.   Since  the  basis  vectors  were 

-T 
P  ,  P,,  and  P^,  d  =  (a3Q  ,a,,  ^34)  =  (0,1,0).   The  new  basis 


vectors  are  PQ,  Px,  P4,  Pa3 


We  find  that  x  3  =  2  >  0  so  we 


go  to  step  8  and  recompute  z.-c,  for  j  =  0,l,.„.,n  . 

The  next  sequence  of  steps  is  lb,  3a,  and  4a,  which 


leads  to  a  new  set  of  s.'s. 

3 


p 

B 

po 

pl 

P2 

P3  p4 

P5 

Pa0 

p  1 

al 

P   O    P   "5 

a2    a3 

po 

bQ-12 

1 

0 

-3/2   1   0 

3/2 

1 

-1 

-1/2   0 

pl 

6 

0 

1 

3/4   3   0 

1/4 

0 

0 

-1/4   0 

P4 

6 

0 

0 

7/4  -3   1 

-3/4 

0 

1 

3/4   0 

Pa3 

2 

0 

0 

l/4©0 

-1/4 

0 

0 

1/4   1 

z.-c. 
3   3 

-2 

0 

0 

-1/4 

^T]  0 

1/4 





s  . 

3 

12 

0 

0 

3/2   2   0 

1/2 

s  . 
3 

8 

0 

0 

10   0 

1 

=  2 


=  -2/-1  =  2 


Steps  la  and  2  bring  P^  into  the  basis  with  the  elimina- 
tion of  P  0  from  the  basis.   The  new  tableau  is: 
a3 


p 

B 

P0 

pl 

P2 

P3 

P4 

P5 

Pa0 

Pal 

P  o 

a2 

P   -3 

a3 

po 

bQ-14 

1 

0 

-7/4 

0 

0 

7/4 

1 

-1 

-3/4 

-1 

pl 

0 

0 

1 

0 

0 

0 

1 

0 

0 

-  1 

-3 

P4 

12 

0 

0 

5/2 

0 

i 

-3/2 

0 

1 

3/2 

3 

P3 

2 

0 

0 

1/4 

1 

0 

-1/4 

0 

0 

1/4 

1 

z.-c. 

3   3 

0 

0 

0 

0 

0 

0 

0 









s  . 

8 

0 

0 

1 

0 

0 

1 
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This  tableau  is  the  final  tableau  since  the  sequence  of 
steps  lb,  3b,  5b,  6a,  and  9b  inform  us  that  we  have  found  an 
optimal  solution  to  the  original  linear  program.   Note  that 
the  complementary  slackness  condition  has  been  maintained, 

n 

that  is,  s«xn  +   y   s.x.  =  0.   Therefore,  the  optimal  solution 
0  0    j^!  3    3 


is 


xQ  =  bQ  -  14  >  0, 


xl  =  X2  =  X5  =  °' 


x-,  =  2,  and 


x4  =  12, 


with  the  optimal  cost  z  =  sn  =  8. 


B 
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V.   PROGRAMMING  TECHNIQUE 
The  linear  programming  technique  described  in  this  thesis 
was  programmed  in  FORTRAN  IV  for  use  on  the  IBM  360/67  com- 
puter.  One  subroutine  is  used  for  the  primal  simplex  itera- 
tion and  another  is  used  for  the  dual  iteration.   The  final 
subroutine  is  used  for  the  addition  of  constraints.   The  main 
(driving)  routine  is  used  to  solve  both  the  full  linear  pro- 
gram and  the  linear  program  using  addition  of  constraints  as 
described  in  this  thesis.   Using  the  main  routine  to  solve 
both  problems  eliminates  any  time  differences  due  to  differ- 
ences in  programming  techniques.   Both  of  the  solution  proce- 
dures were  timed*,  and  the  number  of  iterations  of  each  were 
counted.   Read  and  print  times  were  not  included  in  the 
timing. 


*The  timing  routine  was  developed  by  Lt.  E.A.  Singer,  a 
student  at  the  Naval  Postgraduate  School. 
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VI.   EFFICIENCY  OF  THE  ALGORITHM 
To  eliminate  the  considerable  time  and  effort  required 
to  input  data  by  hand,  a  subroutine  was  designed  which  gen- 
erates random  problems  of  a  large  size.   This  routine  uses  a 
random  number  generator  to  produce  the  elements  of  the  A,  B, 
and  C  matrices.   The  following  criteria  were  used  in  order 
to  insure  the  existence  of  a  bounded  feasible  solution: 

maximize 

n 
z  =   J   c  .x . 

jil   3  ^ 

subject  to 

n 

I 

and 


y   a..x.  -  x  .  =  b.,  i  =  l,...,m, 
S   i]  ]     si     l ' 


c.  <  0.  x.  >  0.  x    >  0,  a..  >  0,  b.  >  0 

:-       J     -  SI-       1J  -       1- 

for  i  =  1 , . . . ,m,  j  =  1 , .  . . ,n , 

l_T_ 

where  x  .  is  the  slack  variable  for  the  i    constraint, 
si 

For  this  investigation  the  subroutine  generated  problems 
having  70  variables  (including  20  slack  variables) ,  20  con- 
straint equations,  and  50  cost  coefficients  using  the  follow- 
ing uniform  distributions: 

a .  .  ,  uniform  (0,1)  ; 

b .  ,   uniform  (0,5)  ; 

c.,   uniform  (-1,0). 

A  total  of  40  problems  were  solved  using  the  random 
problem  generator  described  above.   The  execution  times  for 
these  problems  are  given  in  Table  I  at  the  end  of  this  section 
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The  comparison  of  solution  times  shows  that  all  40  prob- 
lems ran  faster  using  the  method  described  in  this  thesis, 
than  with  the  standard  primal-dual  algorithm.   The  time  dif- 
ferences range  from  15,17  seconds  to  4  8.18  seconds  with  an 
average  time  difference  of  27  e 53  seconds. 
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TABLE 

I 

Addition 

of 

Prob. 

Full 

Array 

Constraints 

No. 
1 

Time  (x.  ) 
(sec. ) 

Iter. 

Time  (y . ) 
(sec.)1 

Iter 

x.  -  y . 

l    7i 

1 

87.82 

23 

56.50 

21 

31.32 

2 

87.84 

23 

56.91 

23 

30.93 

3 

80.22 

21 

56.53 

21 

23.69 

4 

87.93 

23 

56.70 

22 

31.23 

5 

84.06 

22 

58.61 

23 

25.45 

6 

91.63 

24 

56.62 

22 

35.01 

7 

95.45 

25 

58.65 

23 

36.80 

8 

80.31 

21 

56.52 

21 

23.79 

9 

99.60 

26 

65.67 

32 

33.93 

10 

80.16 

21 

56.51 

21 

23.65 

11 

80.18 

21 

56.50 

21 

23.68 

12 

91.65 

24 

65.12 

29 

26.53 

13 

80.19 

21 

57.99 

22 

22.20 

14 

91.74 

24 

64.35 

30 

27.39 

15 

80.20 

21 

56.48 

21 

23.72 

16 

84.02 

22 

57.10 

22 

26.92 

17 

95.52 

25 

59.36 

24 

36.16 

18 

83.95 

22 

58.71 

22 

25.24 

19 

84.14 

22 

65.47 

25 

18.67 

20 

99.25 

26 

70.03 

28 

29.22 

21 

107.08 

28 

58.90 

25 

48.18 

22 

87.85 

23 

56.51 

21 

31.34 

23 

84.01 

22 

61.39 

24 

22.62 

24 

80.16 

21 

56.52 

21 

23.64 

25 

80.17 

21 

56.49 

21 

23.68 

26 

80.16 

21 

56.50 

21 

23.66 

27 

80.19 

21 

56.52 

21 

23.67 

28 

84.21 

22 

58.84 

23 

25.37 

29 

103.27 

27 

57.85 

24 

45.42 

30 

87.83 

23 

59.47 

23 

28.36 

31 

84.07 

22 

59.33 

22 

24.74 

32 

88.19 

23 

56.92 

23 

31.27 

33 

80.28 

21 

56.56 

21 

23.72 

34 

87.96 

23 

56.56 

21 

31.40 

35 

84.11 

22 

57.09 

23 

27.02 

36 

84.05 

22 

56.56 

21 

27.49 

37 

80.20 

21 

56.51 

21 

23.69 

38 

80.18 

21 

56.54 

21 

23.64 

39 

84.03 

22 

68.86 

26 

15.17 

40 

80.29 

21 

56.54 

21 

23.75 

Total 

3454.01 

2352.77 

1101.24 

Average 

86.35 

58.82 

27.53 
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VII.   SUMMARY  AND  CONCLUSIONS 
A  modification  of  the  primal-dual  algorithm  has  been 
presented.   This  modification  difrers  from  the  standard 
primal-dual  algorithm  in  that,  the  constraint  equations  are 
introduced  one  at  a  time,  and  each  subproblem  is  solved  be- 
fore the  next  constraint  is  added. 

This  algorithm  was  programmed  in  FORTRAN  IV  for  the  IBM 
360/67.   The  program  was  designed  so  that  any  given  linear 
program  is  solved  first  by  the  standard  primal-dual  algor- 
ithm, and  then  is  resolved  using  the  modified  primal-dual 
procedure,   Further,  the  same  subroutines  are  used  for  both 
methods  in  order  to  eliminate  timing  bias  due  to  coding 
differences.   In  fact  the  modified  procedure  contains  steps 
which  are  not  included  in  the  timing  of  the  standard  routine 
In  all  cases  the  modified  procedure  was  faster  than  the 
standard  procedure.   Of  4  0  test  problems  the  standard  method 
averaged  86.35  seconds  per  problem,  whereas  the  modified 
method  averaged  58.82  seconds  per  problem.   One  should  not 
judge  the  actual  running  time  of  the  test  problems  since  no 
attempt  was  made  to  improve  the  efficiency  of  the  computer 
program  on  an  absolute  basis;  only  the  relative  speeds  of 
the  two  methods  is  of  importance. 


23 


BIBLIOGRAPHY 

1.   Hadley,  G. ,  Linear  Programming.  Reading,  Mass:  Addison- 
Wesley  Publishing  Co.,  Inc.,  1962. 


24 


APPENDIX  A 


FLOW  DIAGRAMS  OF  THE  COMPUTER  PROGRAM 
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MAIN  PROGRAM 


F0*    of*  wirtf 


P#6£ 


I 


A'g 
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g        />$£§  2  7 
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0- 
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SIMPLEX   ITERATION   SUBROUTINE 
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fa 
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■7 
J   * 


DUAL  ITERATION  SUBROUTINE 
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4.   ADDITION  OF  NEW  RESTRAINT  SUBROUTINE 
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n 
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APPENDIX  B 


FORTRAN  LISTING  OF  THE  COMPUTER  PROGRAM 
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IV  G  LEVEL  1,  MOD  1 


MAIN 


PATC  =  68159 


05/22/31 


00C2 
0003 


0004 

cocs 

0006 
0007 
00C8 
POOO 
0010 
0011 


CO  12 
0013 


UClt 
CO  15 
0016 

0017 
0018 
001Q 
0020 
0021 
0C22 
002? 
0C2<* 
0025 
0026 
0027 
0028 
002° 
0030 
0031 
0032 
0033 
003<* 
0035 
0036 
0037 
0038 
00?o 
OC^O 

0041 

00^2 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
C050 
0C51 
C052 
0053 


0054 
0055 
0056 
0057 

0058 
0059 
0060 


COMMON/ I  NT/ I UNR, IM,M,N, IRAS, ITER 

COMMGN/FLOAT/P,H,C, RT  , ONF,TFM, PINV, 7CR , 7 C , WSR , WS , T , X , ORAR , FP 

NOTE  -  WS(J)  IN  THIS  PROGRAM  AND  S(J)  IN  THESIS  ARE  EQUIVALENT 


EP  =  10.0E-6 
XF  =  URN<0) 
IPROB  =  0 
6  IF( IPR0B.GE.50)  GO  TO  500 
IPROR  =  IPRPR  +1 
CALL  PSEiinn 

M  =  M  «-  1 
N  =  N  +  1 

K0=1  IS  FOR  USF  OF  COMPLETE  TABLEAU 


KO  = 
IM  = 


INITIALIZE  TABLEAU 


9 

10 


11 


zo 


30 


40 


OP 
BT 
BT 
P( 
IT 
DO 
P( 
BT 
3T 
no 

IB 

ON 
TO 
DO 
IF 
oi 

Gr 
PI 

en 

CA 

?c 

zc 

DO 

zr 
zc 
dp 
zc 

DO 

zc 

IN 
WS 
DC 
IF 
WS 
IN 
CO 
wS 
WS 
DO 
WS 


1=2,^ 
1)  =  O.C 
,2)  =  B(  I 


It  1)  =  0.0 


I) 


=  1 
J=l 


N 
1,J»  =  1.0 

(1.1)  =  1.0 

( 1.2)  =  O.C 
3  1  =  1, M 

AS( I )  =  N+I 
E(  I)  =  -1.0 

10  [  =  1 

10  J  =  l 
( I .EO.J) 
NV( I , J) 

TO  10 
NV( I , J) 
NTINUE 
LL  TI"EIT(C, CD 


go 

=  c. 


TO 
0 


=  1.0 


B(  1) 
3(  2) 

11 
B(  2) 
(  1) 

20 
(J) 

20 
(J) 
0  = 
(1) 

30 


=  -1.0 

=  0.0 
1  =  2,  IM 

=  ZC«<2)  -BT( I ,2) 
=  -1.0 
J  =2  ,N 
=  -1.0 
I  =  2,TM 
=  ZC ( J)  -P(T , J) 


=  0.0 

J  =  2,N 

(CD)  .LE.WSt  1)  )G0  TO 
( 1)  =  C( J) 
D  =  J 
NTINUE 

B(  1)  =  WSl  1) 
B( 2)  =  3.0 

40  J  =  2,N 
(J)  =  WS(  1  >  -C(J) 


30 


DETERMINE  VECTOR  TO  INTRODUCE 


50 
60 
70 


IF(  IND.EO.O.OGO  TO  50 

IF( ZC( I  NO)  .OF. 0.0) GO  TO 

CALL  SIMPLX(IND) 

GO  TO  60 

CALL  SIMPLX(l) 

IF(  IUNB-1 ) 7C,130,190 

TFM3  =  9^9999. 


70 


S8888S3 

0000040 
0000045 


OOOnose 

0000055 
0000060 
0010070 
0000080 
0^00080 
0000090 
0000100 


00002*0 
0000265 


0000770 
0000?«0 
0000290 
0000300 
000^305 
0000310 
00003?0 
0000325 
0000330 
0000340 
0000^50 
0000360 
0000370 
0000330 
00003Q0 
0000400 
0000410 
0000420 
0000430 
0000435 
0000440 
0000450 
0000460 
0000470 
000048Q 
00004*50 
0000500 
0000510 
0000520 
000053O 
0000540 
0000550 
000056C 
0000570 
0000580 
00005Q0 
00004O0 
0000510 
0000620 
0000630 


000064" 
0000650 
0000660 
0000670 
00006R0 
C0006QO 
POOO-'OO 
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IV  G  LEVEL  1.  MOD  1 


MAIN 


DATE 


68159 


05/22/31 


0061 
006  2 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 


0076 
0077 
007R 


C079 
0080 
0081 
0082 
0083 
0084 
0085 

0086 


0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 


FIND  THETA  FOR  DUAL  =  MI  N( -S  (  J  ) /ZC  (  J  )  )  OVER  J  FOR  ZCIJXO 

IND  =  0 

DO  80  J  =  1,N 

IF(ZC< J). GE. 0.0)G0  TO  80 

DO  75  1=1, JM 

IF(  IBAS(  I  I  .EC  J)  GO  TO  80 
75  CONTINUE 

XE=  -WS( J)/ZC( J» 

IFITfcMP.LE.XEIGO  TO  80 

TEMP  =  XE 

IND  =  J 
80  CONTINUE 

IF(  IND.EO.O)  GO  TO  100 

CALL  DUAL! IND) 

CALL  SIMPLX(IND) 

GO  TO  60 


CHECK  FOR  INFEASIABILITY 

100  IF(ZCB( 1) )20C,110,120 
110  IF(ZCB( 2) .LT.O.O)  GO  TO 
120  IF( IM.GE.M)  GO  TO  190 

ADO  NEXT  RESTRAINT 

130  IV  =  IM  +  1 

CALL  REST 
IUNB  =  0 


200 


190  CALL  TIMEIT(-1,CL) 

WRITE (6, 5000) I  PR OB, CI , ITE 
.COOOFORMATI  10X  ,  •  PROBLEM • , I  4 , • 
116,'   ITERATIONS',//) 
IF(KO.EO.O)GO  TO  6 

K0=0  IS  FOR  USE  OF  ADDITION  OF  RESTRAINTS 


TIME  IS  •♦-6PF15.6,'   SECONDS  WITH', 


IM  =  2 

K0  =  0 

GO  TO  5 
500  STOP 
200  WRITE(6,4010) 

GO  TO  190 
4010  FORMAT! 10X, 'SOLUTION 

END 


INFEASIBLE' //) 


0000710 
0000720 
0000730 


0000740 
0000750 
0000760 
0000770 
00007*0 
0000790 
0000800 
0000810 
0000820 


0000830 
0000840 
0000850 


0000860 
0000R70 
0000880 
0000890 
0000900 
0000910 
0000920 
0000940 
00O09Q0 


0001000 
0001010 
0001020 
0001030 
0001040 
0001050 
0001070 
0001080 
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IV  G  LEVEL  1,  MOO  1  MAIN  DATE  =  68159  05/22/31 


001  MENS  I  ON  P(53,151),3<52),C(151),BT<53,2), I BAS ( 53 ) , ONE ( 53 ) , TE 
1PINV(53,53),ZCB(2) ,  ZC  ( 1511 ,WSB(?) ,WS( 151) , T ( 2 ) , X ( 53 > , OBAP ( 53 ) 


SUBROUTINE  SIMPLX(J)  0001090 

TEM(2),  0001100 
00011  \r> 

COMMON /t NT /IUNB, IM,M,N,IBAS,ITFR  C001120 

CCM«lON/FLOAT/P,BtC,BTf  ONE.TEM.PINV,  ZCB,  ZC.WSPtWS.T,  X,DBA«,EP  00011  2  5 

IUNB  =  0  0001130 

C 

C    FOR    J=l    THE    ONLY    CHANGE     IS    THE    BASIS    ANO    ZCUI'S 

C 

5  IF(J.EO.l)  GC  TO  130  0001140 

C 

C  COMPUTE  X(J)  ANO  FIND  THETA  FOR  PRIMAL  =  M I N ( X ( f , B ) /X ( I , J ) )  OVER  J 

C     FOP  X(I,J)>0 

C 

T(l)  =  999999.  0001150 

T(  2)  =  T(U  0001160 

DO  20  1=1, IM  0001170 

l  =  iM-i  +  i  ooonao 

X(L)  =  0.0  0001190 

00  10  K=l , IM  0001200 

Y  =  PINV(L ,K)*P(K, J)  0001205 
X(L)  =  X(L)  ♦  Y  0001210 

10  CALL  R0UN0(X(L > ,Y,EP)  0001211 
IF( X(L) .LE.C. )G0  TO  2C  000121^ 
TEM(l)  =  BT(L,1)/X(L)  00012">n 
TEM(2)  =  8T(L,2)/X(L)  0001230 
IFITI  1  1-TEMI 1)  )2C, 11,12  0001240 

11  IF(T(2).LE.TfcM(?) >G0  TO  20  00012^0 

12  T<  1)  =  TEM( 1)  0001260 
T(2)  =  TEMI2)  00012^0 
ID  =  L  0001280 

20  CONTINUE  0001290 

C 

C  COMPUTE  NEW  TABLEAU 

C 

IF(T(  1).FQ. 999999. )G0  TO  50  0001300 

IBAS(IO)  =  J  0001310 

GNE( 10)  =  0.0  0001320 

00  30  1=1 , IM  0001 330 

IF(  I  .EQ.ID)  GO  TO  '0  0001 333 

Y  =  T< 1 )*X( I )  000133? 
BT( I  ,  1)  =  BT( I ,1 )  -  Y  0001334 
CALL  ROUNO (BT( I , 1) ,Y,EP)  0001335 

Y  =  T<  2)*X(  I  )  0001 336 
BT(  I  ,2)  =  BT(I ,2)  -  Y  0001336 
CALL  P0UN0(BT( 1,2) ,Y,EP)  0001337 
XE  =  X( I )/X( 10)  0001360 
DO  30  K=1,IM  0001370 

Y  =  PINV( ID,K)*XE  0001371 
PINV(I,K)  =  PINV(I,K)  -  Y  0001372 
CALL  R0UN0(PIMV( I ,K) ,Y,EP)  000137S 

30  CCNTINUF  00013"5 
BT(ID,1)  =  PT( I0,1)/X( 10)  0001387 
BT(I0,2)  =  BT( 10, 2) /X( 10)  C0013S8 
DO  31  1  =  1, IM  C001386 

31  PINV(IO,I)  =  PINV(  ID, I  )  /X( ID)  0001389 
70  ZCB( 1)  =  0.0  0001390 

ZCB(2)  =  0.0  0001410 

00  80  1=1, IM  0001420 

ZCB(i)  =ZCB(1)  «-ONE(  I)*BT(  I  ,1)  0C01430 

80  ZCBI2)  =ZCB(2)  *ONE ( I ) *"T ( I , 2 )  0001440 

ZC(  1)  =  ONE(  1)  C001450 

OP  90  L=2,N  0001460 

Zf(L)  =  0.  0001470 

DO  90  1=1, IM  0001490 

DP  90  K  =  l  , IM  0001S00 

Y  =  ONF(  I )*PINV(  I  ,K  )*P(K,L)  0001505 
ZC(L)  =  ZC(L)  ♦  Y  0001510 

90  CALL  ROUNO(ZC(L)  ,Y,EP)  0001515 
C 
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FORTRAN  IV  G  LEVEL  1,  MOD  1  SIMPLX  DATE  =  68159  05/22/31 

C  CHECK  FOR  AN  ADDITION  ITERATION  WITH  PRIMAL 


0059 

TEMP  =  0.0 

0060 

DO  120  L  =  1,N 

0061 

IF(WS(L).NE.0.0)  GO  TO  120 

0062 

00  110  1  =  1,  IM 

0063 

IF<  IBAS(  I  )  .EO.DGO  TO  120 

0064 

110 

CONTINUE 

0065 

IF(ZC(L).GE.TEMP)  GO  TO  120 

0066 

J  =  L 

0067 

TEMP  =  ZC(L) 

0068 

120 

CONTINUE 

0069 

ITER  =  ITER  ♦  1 

CO  70 

IF(TEMP.NE.O.O)  GO  TO  5 

0071 

200 

RETURN 

0072 

130 

IBAS(l)  =  1 

0073 

ONE( 1)  =  0.0 

0074 

ID  =  1 

0075 

C 

C  SOLI 
C 

50 

GO  TO  70 

JTION  UNBOUNDED 

0076 

IF( IM.GE.Ml  GO  TO  60 

0077 

IUNB  =  1 

0078 

GO  TO  200 

0079 

60 

WPITE(6,4020) 

0080 

IUNB  =  2 

0081 

4020 

FORMAT! 10X, 'SOLUTION  UNBOUNDED 

0082 

GO  TO  200 

0083 

END 

0001520 
0001530 
0001540 
0001550 
0001560 
0001570 
0001580 
0001590 
0001600 
0001610 
0001615 
0001620 
0001630 
0001640 
0001650 
0001655 
0001660 


0001670 
0001680 
0001690 
0001700 
0001710 
//)  0001720 

0001730 
0001740 
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FORTRAN  IV  G  LFVEL  it  MDU  1 


MAIN 


DATE  =  68159 


05/22/31 


0001 
0002 

C  0  Q  3 
00  C  4 


0^0  5 
0006 
00C7 

ooca 

000Q 
OGIC 
0011 
0012 
0013 
0014 
0015 
0016 
0017 


E 

c 

SUBROUTINE  nilAL(J)  0001750 

001  MENS  I  UN  P(5?,151),H('32l,C(151)fRT(53t2),IBAS(5^),nNE(5^)fTFM(2),  0001760 

lPINV(53t53)(ZC8(2),7C(15l),WSR<2>,WS(151>,T<2),X<51),ORAP<53)  0001 7  7^ 

COMMJN/IN'T/IUNB, IMfM.N.IPAS,ITFR  000178^ 

CrMVON/FLOAT/P,BfC,BT ,ONF ,TFMfPlNV,ZCR,ZC,WSBfWS,TfX,nBAP,FP  0001785 

c 

C  CUMPUTF  MFW  WS(J)  WHICH  IS  EQUIVALENT  TO  S  (  J  > 

C 

XF  =  WS(J)/ZC(J)  0001790 

00  10  K=1,N  0001800 

Y  =  ZC(K)*XF  0001^01 
WS(K)  =  WS<K)  -  Y  000190' 

10  CALL  ROUNO(WS(K)  ,Y, EP1  000180^ 

Y  =  ZCrt(l»*XE  0001004 
WSR(l)  =  WSB(l)  -  Y  0001805 
CALL  R0UN0( WSO( 1 ) ,Y,EP)  C001RQ6 

Y  =  ZCB(2)*XF  0001807 
WSB( 21  =  WSB(2)  -  Y  0001808 
CALL  POUNDIWSTC  >  ,Y,EP)  0001810 
RFTDRN  0001840 
END  0001850 
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FORTRAN  IV  G  LEVEL  1,  MOD  1  MAIN  DATE  =  6*159  05/22/31 


E 


C 

0001  SUBROUTINE  REST  0001860 

0002  ODIfENSION  P( 53, 151 ) , B( 52) , C( 15 1 ) , BT ( 53 , 2) , I RAS ( 53  I , ONE ( 53 ) , TEM( 2 ) ,  0001R70 
1PINV(53,53),ZCB<2) , ZC( 151) , WSB ( 2  I , WS< 15 1 ) ,T< 2 ) , X ( 53 ) , OBAR < 53 )  0001880 

0003  COMMON/INT/IUNB, IM,M,N,IBAS, ITER  0001R°0 

0004  C0MM0N/FL0AT/P,B,C,RT, ONE , TEM, PI NV  ,  ZCB, ZC , WSB, WS , T , X , OBAR , EP  0001  895 
C 

C  COMPUTE  D-BAR  AND  B  VECTOR 
C 

0005  K  =  IM-1  0001900 

0006  BT(IM,1)  =  0.0  0001910 

0007  BT(IM,2)  =  B(K)  0001920 

0008  00  10  1=1, K  0001930 

0009  IFUBASU)  .GT.N)  GO  TO  5  0001935 

0010  L  =  IBASU  )  0001940 

0011  DBAR(I)  =  -P(IM,L)  0001950 

0012  BT(IM,1)  =  BT(IM,1»  ♦  OBAR ( I » *BT ( I , 1 )  0001960 

0013  BT(IM,2)  =  BT(IM,2»  *■  OBAR ( I )*BT ( I , 2 )  0001970 

0014  GO  TO  10  0001<J71 

0015  5  DBAR(I)  =  O.C  0001972 

0016  10  CONTINUE  0001973 

0017  IF(BT( IM,1))5C,2C,30  000198C 

0018  20  IF(BT( IM,2>.LT.0.0)  GO  TO  50  0001990 
C 

C  COMPUTE  NEW  INVERSE  ANO  COST  VECTOR 
C 

30  00  40  J=1,K  0002000 

DO  40  1=1, K  0002010 

40  PINV(IM,J)  =  PINV(IM,J)  «-  OBAR(  I  |*PINV(  I  ,  J)  0002020 

ZC(1)  =  ONF(l)  0002021 

00  45  J=2,N  0002022 

ZC(J>  =0.0  0002023 

00  45  1=1, IM  0002024 

DO  45  IK=1,IM  0002025 

Y  =  ONE(  H*PINV(  I,  IK)*P(  IK,  J)  0002026 

ZC(J>  =  ZC(J)  +  Y  0002027 

45  CALL  ROUND(ZCU)  ,Y,EP)  0002028 

ZCB(l)  =  ZCBUI  -  BT(IM,1)  0002029 

ZCB(2)  =  ZCBI2)  -  PT(IM,2)  0002030 

100  RFTURN  0002031 

C 

C  USED  TO  INSURE  FEASIBILITY  BY  MAKING  THE  R  COMPONENT  NON-NEGATIVE 
C 

50  BT(IM,1)  =  -PT(IM,1)  0002040 

BT(IM,2»  =  -BT(IM,2)  0002050 

DO  60  1=1, K  0002060 

60  DBAR(I)  =  -DPAR(I)  0002070 

PINV( IM.IM)  =  -1.0  0002130 

GO  TO  30  0002140 

END  0002150 
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ORTRAN  IV  G  LtVFL  1,  MOD  1 


MA  IN 


DATF 


6A1S9 


05/22/31 


00  CI 

0002 
0003 
0004 
00C5 
C0C6 


SUBROUTINE  RCUNOU,R,EP) 
A^STARTING  VALUE,  R=ADDFD  VALUE,  EP  =  LOWFST 
IF< A.EU.O.O.rR.^.EO.0.0)  GO  TO  200 
IF(AriS( A/B J.GT.EPI  GO  TO  200 
4  =  0.0 
200  RETURN 
END 


ROUND-OFF  DFSIPFO 


0002160 
0002170 
0002130 
0002190 
0002200 
000221O 
no0?2?0 
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FORTRAN  IV  G  LEVEL  I,  MOD  1  MAIN  DATE  =  68159  05/22/31 


E 


c 

0001  SUBROUTINE  TIME  I T ( N,TI ME )  0002230 
C 

C  N=0  STARTS  CLOCK,  N=-l  STOPS  CLOCK 
C 

0002  IT=N+2  0002270 

0003  GO  TO  (20, 10), IT  0002280 

0004  10  CALL  TIMON(M)  00022<50 

0005  TIMEM=M  0002300 

0006  RETURN  0002310 

0007  20  CALL  TIMOFF(M)  0002320 

0008  TIME=M  0002330 

0009  TIME=(TIMEM-TIME)*26.0  0002340 

0010  RETURN  0002350 

0011  END  0002360 
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FORTRAN  IV  G  LEVEL  1,  MOD  1 


MAIN 


DATE  =  68159 


05/22/31 


00C1 
0002 

0003 
C004 


00C5 
0006 
0007 
0008 
000° 
001C 
0011 
C012 
0013 
OOK 
0015 
0C16 
0017 
0018 
0019 
0C2C 
C021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
C033 


SUBROUTINE  PSEUDC 
ODIMENSION  P<53,15i),fM52),C(151»,BT(53,2),IBAS(53)  , ONE < 53 > , TFM(  2), 
lPINV(53,53),7CB<2),7r<151),WSP<2>,WS(151),T(2),X(5^) , DRAR(53) 

COMMON /I  NT/ I UNR, IM,M,N. IRAS,  ITER 

COMMJN/FLOAT/P,B,C,BT,nNE  ,TEM,PINV,7CB,7C.  ,WSB,WS,T,X,0B4R,EP 

M  =  NUMBER  OF  CONSTRAINTS,  N  =  NUMBER  OF  VARIABLES  INCLUDING  M  SLACK 

VARIABLES 
MAXIMUM  M  =  50,  MAXIMUM  N  =  150 

CHANGE  STATEMENTS  1  AND  2  TO  DESIRED  PROBLEM  SI ZF 


10 


20 


30 


40 


50 


M  =  20 
N  =  70 


IM  = 
IN  = 
Kl  = 
K2  = 
C  (  J  » 


♦  1 

♦  1 

-  M  + 
♦  I 


1 


GO  Tn  10 


=  -URN( 1) 
no  20  J=2,K 
IF(C( J) .LE.-.001 ) 
C( J)  =  0.0 
00  20  1=2, IM 
P(  I,  J)  =  URM1) 
IF(P<  I  ,J).GT. 0.001 >  GO 
P(  I  ,J)  =  0.0 
CONTINUE 
10  30  1  =  1  ,  M 
H(  I  )  =  URN(  1)*5.0 
IE(R( I ) .GT.0.C1)  GO  TO  30 
R( I )  =  0.0 
CONTINUE 
DO  40  J=K2,IN 
C(J)  =0.0 
DO  40  I  =  2, IM 
P( I , J)  =  O.C 
DO  50  1=2,  IK 
K  =  N-  M  +  I 
P(  I  ,K)  =  -l.C 
RETURN 
END 


TO  20 


C00237Q 
0002380 
0002290 
0002400 
D002410 


0002430 


0002450 

0002470 
0002480 
0002490 
0002500 

0002520 
OO32530 
0002^40 
0002550 


0002580 


0002640 
0002650 
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